Phase-space methods in nuclear reactions around the Fermi energy 
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Some prescriptions for in- medium complex particle production in nuclear reactions are proposed. 
They have been implemented in two models to simulate nucleon-nucleus (nIPSE) and nucleus- 
nucleus (HIPSE 1 ) reactions around the Fermi energy 0, Our work emphasizes the effect of 
randomness in cluster formation, the importance of the nucleonic Fermi motion as well as the role 
of conservation laws. The key role of the phase-space exploration before and after secondary decay 
is underlined. This is illustrated in the case of two debated issues: the memory loss of the entrance 
channel in central collisions and the (iV, Z) partitions after the pre-equilibrium stage. 
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In nuclear reactions around the Fermi energy, nuclei 
can break into several pieces of various sizes: the so- 
called multifragmentation process |3|. A striking experi- 
mental feature is the large number of accessible charge 
and energy partitions. In order to understand the statis- 
tical aspects of the explored phase-space, several scenarii 
have been proposed. However, in view of the complexity 
of nuclear reactions due for instance to impact parameter 
mixing, pre-equilibrium emission and thermal decay, it is 
hard to trace back the process of cluster formation. This 
issue remains a highly debated question. During the past 
three years, we have try to answer the following question: 
"What are the minimal hypothesis required to reproduce 
the characteristics of nuclear fragmentation?" . Based on 
a detailed and a successfull comparison with experimen- 
tal data, we have ended up with surprising conclusions. 
Simple prescriptions have been found for the formation 
and the emission of complex particles 0, 0, H- First, 
these rules are introduced and illustrated. In a second 
part, we discuss two aspects of current interest related to 
the phase-space exploration after pre-equilibrium stage. 



I. RULES FOR THE FORMATION AND THE 
EMISSION OF CLUSTERS 
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1 The computer program of the model is available at 
http://caeinfo.in2p3.fr/theorie/theory_lacroix.html 



FIG. 1: Schematic representation of a three step nucleon- 
induced reaction. A nucleon with energy close to the Fermi 
energy is first absorbed by a nucleus. Then two steps are 
identified for pre-equilibrium emission: the formation of the 
cluster and its emission. 

In this section, rules for the cluster formation and 
the production of fragmentation partitions are defined. 
These rules have been extensively discussed in ref. Q 
and we only give here a short summary. Let us con- 
sider a proton colliding a heavy target with an energy 
Eb close to the Fermi energy. A schematic three steps 
scenario is considered (see figure [Q. First, the incident 
nucleon is absorbed. The second step corresponds to the 
in- medium formation of the cluster. We assume that 
clusters are made of nucleons randomly picked inside a 
Thomas-Fermi distribution. Cluster properties, position 
r c and the kinetic energy Sk are calculated on the basis of 
the characteristics of the nucleons in phase-space. Using 
a Monte-Carlo sampling, an ensemble of configurations 
(called "total accessible phase-space") is obtained. As an 
illustrative example, the phase-space of an a particle is 
given in top right of figure [3 The last step corresponds 
to cluster emission. All configurations obtained in the 
second step will not necessarily lead to the emission of a 
cluster. Two constraints on the emission can be identi- 
fied. The first one, which is independent of the entrance 
channel, is due to the mutual interaction between the 
cluster and the other nucleons. Given an interaction po- 
tential between the cluster and the emitter, denoted by 
Va-\-a c (t c ) : the emitted particle must overcome the bar- 
rier, denoted by Vb- This gives a lower limit on the 
cluster kinetic energy. The second constraint is directly 
dependent on the reaction type and is due to the energy 
balance. The constraint of a positive excitation energy 
for the configuration leads to 



Eb — Q — V a +a c (tc) > £k(r c ) 



(1) 



where Q is the Q-value of the configuration. Therefore, 
only a fraction of the total phase-space accessible for the 
cluster will indeed lead to emission in the continuum. 
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cluster size and multiplicity are not fixed a priori, there- 
fore a specific method should be developed to describe 
the nucleosynthesis (i.e the aggregation of nucleons) at 
the early stage of the reaction. Finally, the possible 
secondary decay should be accounted for. A guideline 
for practical implementation of these different aspects is 
given m ref. g. Two models (called n-IPSE 1 2] and 
HIPSE 2 [l| based on the very same assumptions have 
been developped with a very good agreement with a large 
body of experimental data. For instance, in nucleus- 
nucleus reactions, charge distributions, average kinetic 
energies as well as the fluctuations of kinetic energies, 
angular distribution, parallel velocities, event by event 
internal correlations, isospin diffusion effects, have been 
properly reproduced in Xe+Sn reactions for beam energy 
from 25 to 80 MeV/A. In the models, few parameters 
depending only on beam energy have been adjusted by 
comparison with data. These parameters are not changed 
when the target and/or the projectile are changed. This 
aspect is essential to insure the predictivity of the models. 



FIG. 2: Top- left: Two-body potential between the a and 
the emitter. Above the line (I), the cluster cannot be emitted. 
This upper limit is directly given by the energy balance of the 
reaction. Below the line (II), the cluster cannot overcome the 
barrier (since here, quantum tunnelling is not considered). In 
between the two lines, there is a small "escape window" for 
the emission of the cluster. Top-right: Total available phase- 
space of the cluster. This latter is significantly reduced due to 
the energy constraint. The two curves correspond respectively 
to the lower and upper limit in the kinetic energy. Bottom: 
Calculated kinetic energy distribution (open squares) of the 
a particle obtained by propagating each configuration in the 
"escape window" up to infinity. The calculated spectrum is 
compared with the experimental data (black circles). 



This fraction corresponds to the "explored phase-space". 
These two constraints are shown in Fig. [2| (top left) for 
a proton-induced reaction at Eb = 39 MeV. There, an a 
particle can only be emitted in a small interval of kinetic 
energy (called "escape window" in the following) leading 
to a significant restriction in phase-space. According to 
the energy constraint, all configurations between the two 
lines lead to the emission of an a particle (Figure (top- 
right)). 

If we now apply the rules and propagate classically all 
explored configuration with the interaction potential, we 
obtain a distribution of emitted kinetic energy directly 
comparable to experiments (bottom of figure EJ). It is 
observed that a direct application of the rules already 
leads to a very good agreement with data [5|. 

However the simplified three steps scenario described 
above allows only for qualitative comparisons with data. 
In order to provide quantitative comparisons, the follow- 
ing effects should be incorporated: the influence of the 
impact parameter and the related geometric effects, the 
in-medium nucleon-nucleon collisions which distort the 
Fermi motion of nucleons. It is worth to notice that the 



II. PHASE-SPACE EXPLORATION: TWO 
ILLUSTRATIONS 

A detailed comparison of the models with the experi- 
mental data can be found in [l| for nucleus-nucleus and 
for nucleon-nucleus reactions. Besides the fact that 
comparisons are made on the phase-space after the sec- 
ondary decay, HIPSE and nIPSE models give also access 
to configurations before secondary emission. This inter- 
mediate stage is of particular interest for making the link 
with statistical physics. In the following, we give two ex- 
amples where the models provide important information 
in the nuclear context. 



A. Loss of memory of the entrance channel in 
central events 

A possible way to trace back the complete loss of mem- 
ory of the entrance channel is to consider the degree of 
shape equilibration. Following ref. p, we define the 
isotropic ratio as (Ri so ) = E±/2E// where E± and Ej / 
correspond respectively to the average total transverse 
and parallel energy. For a spherical symmetric partition, 
we do expect (Ri SO ) = 1, this limit corresponding to a 
complete shape equilibration. In figure (0), we have plot- 
ted the evolution of (Ri SO } as a function of multiplicity 
for the reaction Xe+Sn at 50 MeV/A. We observe an in- 
crease of this quantity and a saturation around a value of 
0.75. A similar behavior has been found in experiments 



1 n-IPSE: nucleon-Ion Phase-Space Exploration 

2 HIPSE: Heavy-Ion Phase-Space Exploration 
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and assigned as a incomplete loss of the entrance chan- 
nel memory. In that case, the multiplicity has been as- 
sumed to be strongly correlated with the centrality of the 
reaction. In HIPSE, the average multiplicity indeed in- 
creases with the centrality. However, large fluctuations 
around the mean value are also observed. As a conse- 
quence, a given multiplicity corresponds to a large set of 
impact parameters (denoted by b). The great advantage 
of HIPSE is that it gives direct access to b. In the left 
side of figure ©, the evolution of (Ri so ) is presented as 
a function of b. We clearly observed a saturation around 
1 for smallest impact parameters. This observation leads 
to two conclusions for the model. First, the multiplicity is 
not a perfect selector of the centrality and the saturation 
below 1 is a spurious saturation due to impact param- 
eter mixing. Second, most central events in HIPSE are 
completely equilibrated in shapes and a complete loss of 
memory of the entrance channel is observed. 



1.25 




FIG. 3: Evolution of the isotropic ratio as a function of the 
total multiplicity (left) and impact parameter (right). The 
black circles indicate the average evolution. 



B. Isospin equilibration 

With new radioactive beams facilities, we do expect 
to explore a large set of nuclei in the nuclear chart. In 
HIPSE, a particular care has been taken to treat consis- 
tently the isospin degrees of freedom. For the moment, 
HIPSE has only been applied to reactions with stable 
nuclei but can similarly be used for radioactive beams. 
This aspect is illustrated in the reaction 48 Ca+ 40 Ca at 
beam energy E B = 25 MeV/A. In figure (@J), the (TV, Z) 
distributions before the decay respectively for the quasi- 
projectile (QP) and the quasi-target (QT) in the left side 
and for the mid-rapidity clusters and compound nucleus 
in the right side are shown. This figure clearly shows the 
partial isospin equilibration between the QP and QT, 
while clusters in the mid-rapidity region or close to com- 
pound nucleus are equilibrated in average. It also illus- 
trates the large fluctuations in the nuclear chart. 



48Ca+40Ca - E B =25 MeV 




FIG. 4: Illustration of the Neutron-Proton phase-space explored for the Quasi- Target (QT), Quasi-Projectile (QP) in the left 
side and for the Mid-rapidity clusters and Compound nucleus in the right side for the reaction 48 Ca+ 40 Ca at beam energy 
Eb — 25 MeV/ A (all impact parameters are mixed). In both cases, the largest area indicates the total phase-space available 
(in the mass table used in HIPSE), points indicates nuclei in the valley of stability, dashed lines correspond to the Target and 
Projectile N/Z and the central solid line to the equilibrium N/Z. 



III. CONCLUSION 



We propose simple prescriptions that may be used to 
describe the pre-equilibrium emission of clusters in the 
course of nuclear reactions. These rules are based on 
a random sampling of the nucleons taking into account 
the Fermi motion and a proper account of nuclear ef- 
fects as well as the conservation laws. Using these rules, 



two models dedicated respectively to nucleon nucleus 
(nIPSE) and nucleus-nucleus reactions (HIPSE) around 
the Fermi energy have been developed 

Besides the good reproduction of a large variety of ex- 
perimental observations on a large set of different reac- 
tions and energies, these models underlines the impor- 
tance of phase-space explored before and after the sec- 
ondary decay. This aspect is illustrated in two exam- 
ples of present interest, namely the memory loss in most 
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central collisions and the (TV, Z) partitions after the pre- equilibrium stage. 
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